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ABSTRACT 

We present a determination of the 'Cosmic Optical Spectrum' of the Universe, i.e. the ensemble 
emission from galaxies, as determined from the red-selected Sloan Digital Sky Survey main galaxy sample 
and compare with previous results of the blue-selected 2dF Galaxy Redshift Survey. Broadly we find 
good agreement in both the spectrum and the derived star- formation histories. If we use a power-law 
star-formation history model where star-formation rate esc (l+z)^ out to z = 1, then we find that of 2 to 
3 is still the most likely model and there is no evidence for current surveys missing large amounts of star 
formation at high redshift. In particular 'Fossil Cosmology' of the local universe gives measures of star- 
formation history which are consistent with direct observations at high redshift. Using the photometry 
of SDSS we are able to derive the cosmic spectrum in absolute units (i.e. W Mpc^'^) at 2-5A 
resolution and find good agreement with published broad-band luminosity densities. For a Salpeter IMF 
the best fit stellar mass/light ratio is 3.7-7.5 Mq/Lq in the r-band (corresponding to ^Istarsh = 0.0025- 
0.0055) and from both the stellar emission history and the Ha luminosity density independently we find 
a cosmological star-formation rate of 0.03-0.04 h Mq yr~^ Mpc"'^ today. 

Subject headings: cosmology: miscellaneous, observations - stars: formation 



1. INTRODUCTION 

The comoving star-formation rate (SFR) of the Universe 
is in decline. Since z = 1 it has dropped by a factor of 3-15 
(Cowie et al. 1999; Lilly et al. 1996). At higher redshifts it 
may have been constant or declining from z = 1 to z = 5 
but the evidence is that z = 1 is a critical point for the 
universal average SFR (Madau et al. 1996; Steidel et al. 
1999). 

This is a very dramatic conclusion and worth tackling 
with a variety of complementary observational techniques. 
One must of course consider how SFR is measured at 
different redshifts: for z > 1 the principal observational 
probe has been rest-frame ultraviolet emission since it is 
easily accessible in optical bandpasses. Because the UV 
derives from young stellar populations (t << IGyr for 
Arest < 2000A) the UV flux produced per galaxy will be 
proportional to the SFR. There are details to do with the 
contamination from older populations but for z > 1 and 
Arest < 2000A these arc minor (Madau et al. 1998). A 
very important issue though is dust extinction which can 
be many magnitudes in the UV and in fact represents the 



principal uncertainty in whether the SFR drops off again 
at high redshift (see for example Pettini et al. 1998; Steidel 
et al. 1999). A comprehensive discussion of UV dust ex- 
tinction issues is given by Bell (2002). Another important 
issue is the treatment of cosmological surface-brightness 
dimming which may result in surveys missing a dramatic 
increase in SFR at high redshift (Lanzetta et al. 2002). 

At lower redshifts, where the rest-frame optical spec- 
trum is accessible, alternative methods of calculating SFRs 
from line emission can be used. The simplest is Ha which 
traces the number of Lyman continuum photons (Glaze- 
brook et al. 1999; Hopkins et al. 2000). This technique 
has the advantage of using radiation emitted at red wave- 
lengths and is thus considerably less sensitive to dust ex- 
tinction compared to the UV. Other lines such as H/3 
(Tresse & Maddox 1998a) and [OH] (Hogg et al. 1998) are 
used but their relationship to SFR is more complicated. 
For z < 1.5 some of these lines are accessible in the opti- 
cal and near-infrared windows and can be used to measure 
SFRs. A review is given by Hogg (2002) — the vast ma- 
jority of the surveys that have been carried out indicate 
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derived from direet high redshift measurements via lumi- 
nosity densities. 

The Sloan Digital Sky Survey (SDSS, York et al. 2000) 
provides many advantages over the 2dFGRS for this type 
of analysis: it is of higher spectral resolution (though we 
are not yet able to exploit this in this paper) and will be 
about four times larger upon completion. The spectral 
wavelength coverage is larger (3600A < A < 8OOOA for 
2dFGRS and 3800A < A < 9200A for SDSS). The photo- 
metric calibration is much better as each SDSS multi-fiber 
observation contains numerous standard stars and can be 
individually calibrated, whereas for the 2dFGRS a mean 
calibration was applied to all the survey spectra. The ex- 
cellent spectrophotometry is borne out by the good agree- 
ment between synthetic colors computed from the spec- 
tra and actual colors which agree, in average spectra to 
< ±5% (Tremonti 2002, private communication). Also 
the accurate five-color photometry allows comparisons of 
photometric constraints with color constraints. In partic- 
ular SDSS is selected in the r band (A = 6200 ± 600 A), 
whereas 2dFGRS is selected in bj (A = 4600 ±700 A). Thus 
one would expect SDSS to be more biased toward old mas- 
sive galaxies and 2dFGRS to be biased toward yomig, star- 
forming galaxies. The comparison between the two allows 
us to investigate the uncertainties in the determination of 
the SFH. 

In this paper we compute the cosmic spectrum for the 
SDSS local volume, we make a direct comparison with that 
derived from 2dFGRS by BG02, and we derive new con- 
straints on star-formation history models from the SDSS 
spectrum. The plan of this paper is as follows. In Section 2 
we describe the SDSS data and our methods for combining 
the spectra to form cosmic spectra. In Section 3 we de- 
scribe our modeling and fitting procedure and the outcome 
of the comparison of best fitting SFHs between the SDSS 
and 2dFGRS surveys. We also test the consistency of mod- 
els of SFH which form a lot of stars at 2: > 1. In Section 4 
we generate an absolute cosmic spectrum which we show 
in physical units. We use this to estimate emission- line 
luminosity densities and the current SFR density. Finally 
we give our summary and conclusions (Section 5). 

Throughout this paper we take Hq = 70kms~^ Mpc~^, 
^mo = 0-3 and = 0.7 for our cosmological quantities, 
and where appropriate, define h = iIo/100kms~^ Mpc~^. 

2. SAMPLES USED IN THIS ANALYSIS 

The Sloan Digital Sky Survey is a digital CCD survey 
in 5 optical bands which intends to cover up to 10 000 
deg^. An overview is given by York et al. (2000). The 
imaging camera is described by Gunn et al. (1998), the 
ugriz photometric system and calibration by Fukugita 
et al. (1996), Lupton et al. (1999), Hogg et al. (2001) 
and Smith et al. (2002). A large fraction of SDSS data 
is currently available to the entire astronomical commu- 
nity (Stoughton et al. 2002). The coordinate system is 
defined to a precision of better than 0.1 arcsec (Pier et al. 
2002). The main galaxy sample is essentially a magnitude 
limited spectroscopic sample. It is selected as virtually all 
galaxies in the photometric area with a Petrosian mag- 

In reality we can only measure galaxies down to some limiting magnitude, so any calculated cosmic spectrum is just an estimate of the true 
Cosmic Spectrum. 

There is of course AGN activity but this is a negligible contribution to the total optical emission as we will see later 



considerable evolution in the SFR for z < 1 although there 
is disagreement on the amount of evolution required. 

A complementary method has been to probe the far- 
infrared emission of galaxies (10/zm < A < 300/im) where 
dust-processed UV is re-emitted thermally. At high red- 
shifts this must be followed into the sub-mm radio wave- 
lengths. The emission has been used to constrain SFR at 
low-redshift (Rowan- Robinson et al. 1997) and at high red- 
shift (Hughes et al. 1998) but this approach suffers from 
both micertainty in the dust modeling and a lack of spec- 
troscopic redshifts. The latter issue is addressed by fitting 
the spectral energy distributions (SED) with a series of 
templates in order to photometrically estimate the red- 
shifts, however there are huge degeneracies between the 
photometric redshift estimate and the assumed tempera- 
ture of the dust SED (Blain et al. 2002). 

All these approaches involve estimating a luminosity, ei- 
ther continuum or line, per galaxy and then multiplying by 
the space density in order to give a luminosity per comov- 
ing volume. At this point the scale factor, SFR/luminosity, 
which is where the main uncertainties arise, allows trans- 
formation to SFR per unit volume. The light budget per 
volume is a useful quantity because it allows the stellar 
emission history of the Universe to be decoupled, in a 
sense, from its dynamical history (i.e. changes in the num- 
ber of counted objects by processes such as galaxy forma- 
tion and galaxy-galaxy merging). This use of luminosity 
density is a direct method, in which an observed luminosity 
density at a given redshift is converted to a SFR density 
at the same redshift. 

An alternative approach is that of fossil cosmology 
where the past history of the Universe is determined from 
its current contents. This can be done by examining the 
resolved stellar populations in the Local Group (Hopkins 
et al. 2001) or in ensembles of galaxies, for example in 
early type galaxies (e.g. Bernardi et al. (2002); Eisen- 
stein et al. (2002)). Our approach is to look at the en- 
semble of all galaxies; the 'Cosmic Optical Spectrum' of 
the local Universe. This represents the luminosity-scaled 
spectra summed over all galaxies. The cosmic spectrum 
can be thought of as the total emission from all the ob- 
jects in a representative volume of the Universe. Objects 
contribute to the cosmic spectrum according to their lu- 
minosity. As in the case for an individual galaxy, this 
spectrum contains a luminosity-weighted mix of features 
from both old and young stars^^ and we can fit models 
of star- formation history to it. In particular because the 
cosmic spectrum represents an average, it will represent 
the end point of the average SFH. Thus we can fit much 
simpler models to the cosmic spectrum than are required 
for the spectra of individual galaxies, since we expect the 
SFH history of the Universe, as a whole, to vary smoothly 
with time. 

This ensemble approach was applied by Baldry et al. 
(2002, hereafter BG02) to the cosmic spectrum (meaning 
the optical spectrum per unit volume) of 166 000 galax- 
ies in the 2dF Galaxy Redshift Survey (2dFGRS, Col- 
less et al. 2001) and derived constraints on allowable star- 
formation histories (SFH) which agreed well with results 
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nitude r < 17.77. The overall targeting completeness is 
92% (Blanton et al. 2002a). For galaxies in this magni- 
tude range the Petrosian magnitude is close to total, they 
are also close to 'model magnitudes' (which are obtained 
by profile fitting) which represent an alternative method 
of trying to estimate total magnitudes. Ninety-eight per- 
cent of the galaxies span a redshift range of < z < 0.25 
with a median redshift of 0.10. Full details of the spec- 
troscopic main galaxy sample are given by Strauss et al. 
(2002). The sample wo consider in this paper uses a total 
of 153000 galaxies selected to r < 17.70 (magnitudes are 
from the PHOTO v5.2 software corrected for Milky Way 
reddening; see Stoughton et al. (2002) for details) from 
the SDSS spectroscopic data. The sample is data from 
463 survey quality plates in the North Galactic Cap which 
were observed between MJD 51433 and 5235 and 86000 
galaxies lie in the range 0.01 < z < 0.11, the principal 
range considered in this paper. Only galaxies with secure 
redshifts were selected. 

The spectra are taken through 3.0-arcsec diameter fibers 
and have a wavelength range of 3800-9200A and a spec- 
tral resolution A/AA = 1800 which is approximately con- 
stant across the spectrum. The signal-to-noise ratio is 
> 4 per pixel (pixels width ~ 1-2 A). We use 'smear' cor- 
rected spectra reduced using SPECTR02D software v4.9. 
The smear procedure uses dithered telescope pointings to 
make a low-order correction to allow for aperture effects, 
it changes the large scale spectral response at the < ±2% 
level in the mean colors. 

Following BG02, the spectra are combined in redshift 
slices (using the redshifts from the Princeton ID spectra 
pipeline, Schlegel et al. (2003)) by scaling them to their r 
luminosity and summing them. This scaling provides for 
a first-order aperture correction. Thus in each slice we get 
a total luminosity spectrum for all galaxies down to the 
luminosity limit of the survey at that redshift. Because of 
the good spectrophotometry of SDSS it is not necessary 
to apply a correction as was done by BG02. 

There are two important issues in using the cosmic spec- 
trum to constrain star- formation histories, the first is lumi- 
nosity bias and the second is aperture bias. These redshift 
surveys represent the spectra of all galaxies down to a fixed 
apparent magnitude. For low redshift galaxies the mag- 
nitude limit will correspond to a faint luminosity. Since 
luminosity densities from typical galaxy luminosity func- 
tions (approximately Schechter (1976) functions) tend to 
converge for M > M*, where M* is the Schechter 'break' 
absolute magnitude; once we have sampled more than 2-3 
mags below the observed cosmic spectrum would be 
very close to the total spectrum. To be quantitative, for 
a faint-end slope in the range —1.3 < a < —1.0 (Blanton 
et al. 2001; Cole et al. 2001; Cross & Driver 2002; Madg- 
wick et al. 2002), integrating down to M* -|- 3 shows that 
between 84 94% of the light has been sampled. Thus we 
expect samples to be highly complete in luminosity at low 
redshift and give a similar result whether they are selected 
in r or 6j. 

At high redshift we expect the blue/red sample selection 

to become most prominent in the luminosity bias. In par- 
ticular for the redshift range we consider the SDSS red se- 



lection does not sample below the 4000 A break (Figure 1) 
where the UV from young stellar populations dominates. 
Thus we would expect the relative bias to run in favor of 
older SFHs in the SDSS sample as we approach z ~ 0.1 
slices. 

In contrast, at low redshift we expect aperture effects to 
be most apparent since the fixed size in arcseconds of the 
spectroscopic fiber aperture corresponds to smaller phys- 
ical scales in the galaxy. At high redshift the effects are 
reversed the fiber apertures ought to sample most of the 
light of a galaxy, but we are only seeing the most luminous 
galaxies and red versus blue selection becomes more im- 
portant. These sample biases are quantified in Figure 2 
which compares the two surveys. 

The SDSS data gives us an opportunity to test for aper- 
ture bias in a way which was not possible for the 2dFGRS 
work: because of the excellent SDSS imaging and photom- 
etry we can use galaxy colors as a function of aperture size 
as a diagnostic of aperture bias. 

For comparison between the two samples at a common 
redshift, we define the redshift limits so that the surveys 
penetrate to the same relative depth below M*. For ex- 
ample in the low redshift volume we use 0.02 ^ z < 0.05 
for SDSS and 0.03 < z < 0.06 for 2dFGRS, the redshift 
upper limit corresponding to M* -|- 2.6 in both cases. The 
2dFGRS cosmic spectra are taken from those published by 
BG02. In Table 1 we define low, medium and high red- 
shift ranges (A, B and C) for which SDSS and 2dFGRS 
are approximately equivalent in luminosity range. Because 
the 2dFGRS sample is slightly deeper than the SDSS sam- 
ple, this correction also works in the direction of minimiz- 
ing the physical aperture difference, since SDSS fibers are 
bigger (in arcseconds) than 2dFGRS fibers. This proce- 
dure reduces the difference in physical aperture diameter 
in ranges A,B,C from ~50% to ~25%. 

Figure 3 shows a comparison between the normalized 
cosmic spectra for SDSS and 2dFGRS for the volumes 
A,B,C. The 2dFGRS cosmic spectrum has been spec- 
trophotometrically corrected against the model fits as de- 
scribed by BG02. The SDSS spectrum uses the native 
spectrophotometry. Despite the differences in sample se- 
lection and spectral resolution the cosmic spectra look very 
similar, especially in the absorption features which we use 
to constrain the SFH. To be quantitative, in region A 
where we expect the selection bias to be smallest, the RMS 
of the ratio of the two spectra over the wavelength range 
4000A-8000A is 2.6% (comparable to the uncertainties in 
the spectrophotometric modeling quoted by BG02), after 
smoothing the SDSS spectrum to 2dFGRS resolution. We 
can break this down into low-pass (i.e. continuum) and 
high-pass (line) variations using a 200A smoothing filter. 
For the smoothed spectra the RMS ratio is 2.4% whereas 
in the spectra with the smooth component divided out the 
RMS ratio is only 0.7%. They key point is the high-pass 
stellar absorption line information is almost identical in 
the two spectra. The similarity of the cosmic spectra 
from the two surveys indicate the selection difference is 
minimal at low redshift. Comparing with the uncertain- 
ties qiioted by BG02 (which arc mostly systematic) we 
would expect our SFHs from SDSS to agree at the 95% 



Another very low resolution comparison is provided by the CIE (1986) chromaticity values of the spectrum, we get (0.340, 0.340) very close 
to the BG02 value of (0.345,0.345) and still close to white 
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Fig. 1. — Comparison of the selection functions of the two surveys. The SDSS r selection filter (solid lines, Stoughton et al.) and 2dFGRS 
bj selection filter (dashed lines, Hewett & Warren private communication, these include atmosphere and telescope) are shown at redshifts 0.0 
and 0.2 (bracketing the redshift range used in this paper). At high redshifts, the bj filter penetrates below the 4000A breaJc and is therefore 
more sensitive to younger stellar populations. The dotted line is a theoretical cosmic spectrum from our analysis. The bj curve was converted 
to photon-response equivalent for comparison with the r curve. 




Fig. 2. — Comparison of sample bias in the SDSS and 2dFGRS redshift surveys. Left: Luminosity bias versus redshift. This shows 
the fraction of total light of the galaxy population sampled down to the survey magnitude limit, calculated by integrating the appropriate 
Schechter luminosity function for each survey (Blanton et al. 2001; Madgwick et al. 2002, both are consistent with a faint-end slope of —1.2). 
Right: Projected aperture size versus redshift. 

Table 1 

Luminosity-bias equivalent redshift ranges (low/medium/high) between the SDSS & 2dFGRS samples 



Region 


2dFGRS redshift range 


SDSS redshift range 


hmiting magnitudes" 


A 


0.025 <z< 0.06 


0.015 <z< 0.05 


M* -1- 4.5 to -1-2.6 


B 


0.06 <z< 0.10 


0.05 <z< 0.08 


+ 2.6 to +1.4 


C 


0.10 <z<0.13 


0.08 < z < 0.11 


A'h + 1.4 to +0.7 



" Approximate range of limiting absolute magnitudes from the low-redshift cut to 
the high redshift cut relative to the Schechter-break luminosity 
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Fig. 3. — Comparison of SDSS and 2dFGRS cosmic spectra. Three spectra are shown for regions A, B and C as defined in Table 1. The 
lower set (thin lines) are for SDSS and the upper set (thick lines) are for 2dFGRS (spectro-photometrically corrected). All are normalized by 
setting the mean flux between 4200 and 5800A to unity but are offset by multiples of 0.2 for clarity. Note the resolution of SDSS is two to 
three times better than 2dFGRS. 
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confidence level with those of BG02. It is this detailed 
modeling of the SFH, to which we now turn. 

3. STAR FORMATION HISTORIES DERIVED FROM COSMIC 

SPECTRA 

3.1. Methods 

For our modeling of star-formation history we use the 
empirical 'double power-law' parameterization used by 
BG02. The star-formation rate (SFR) is a simple func- 
tion of redshift with a break at redshift unity: 

(1 + z)'^ for z < 1 
SFR a (1 + z)" for 1 < 2 < 5 

for z > 5 

The two power-laws are normalized at z = 1 and star- 
formation is started at z = 5. This simple fitting form is 
chosen because it already provides a good match to the 
range of observations of cosmic SFH (for example, see fig- 
ure 9 of Steidel ct al. 1999), has a small number of parame- 
ters and provides acceptable fits. Other parameterizations 
are possible, see BG02 for some examples. All are princi- 
pally measuring the ratio of old to young stars weighted 
in some fashion. 

The fitting of star-formation histories to the cosmic 
spectra proceeds using standard evolutionary synthesis 
techniques and follows that of BG02 for the 2dFGRS 
data. We use the PEGASE.2 evolutionary synthesis mod- 
els (Fioc & Rocca-Volmerange 1997). We assume a uni- 
versal Initial Mass Function (IMF) which is independent 
of cosmic epoch and for this we use the IMF of Salpeter 
(1955). Evolutionary tracks are from the 'Padova group' 
(Bressan et al. 1993) and the stellar atlas comes princi- 
pally from the theoretical stellar spectra of R. L. Kurucz 
as given by Lejeune et al. (1997). Metallicity evolution in 
this code follows the prescription of Woosley & Weaver 
(1995). The fiducial extinction was chosen as that for 
a inclination-averaged disk geometry but this makes lit- 
tle difference because, as we will see below, our analysis 
primarily uses the high-pass spectral information and we 
varied the extinction to test the effect on the broadband 
color information. 

We use an approach of consistent chemical evolution 
where the interstellar medium for forming new stars is 
continuously enriched by the death of old stars. Using 
the PEGASE.2 closed box model, we control this by a 
SFR normalization r defined by BG02. This ranges from: 
r ~ 0.3 where only a fraction of the available gas is used to 
form stars so that the metallicity remains low, to; r ~ 1.4 
where a total mass of stars is formed over time greater than 
the mass of gas initially available so that chemical evolu- 
tion is significantly faster. The metallicity monotonically 
increases with time and different values of r correspond ap- 
proximately to different end-point metallicities. We quote 
end-point metallicity values from PEGASE.2 which are lu- 
minosity weighted. We implicitly assume that the range 
of metallicities at each galaxy epoch can be represented 
by an average metallicity for the purposes of evolution- 
ary synthesis. This approach is more complicated than 
using simple-stellar populations (constant metallicity) but 
is consistent with the logical assumption that older stellar 
populations should, on average, have a lower metallicity 
than younger populations. 



One approach to fitting observed spectra to models 
would be to compute line indices for both data and mod- 
els (e.g. Kauffman et al. (2002)). The approach we use 
here is to fit the entire spectrum using a statistic. For 
establishing our goodness of fit, we follow BG02 and com- 
pute different figure of merits (FOMs) from high-pass 
and low-pass filtered data. As in BG02 we calculate the 
figure-of-merit quantities by summing over all wavelengths 
except near strong nebular emission lines (Oil 3727A; OIII 
5007A; Ha 6563A; Nil 6583A; SlI 6716A & 6730A) as we 
are only interested in the stellar emission. The first step 
is to fit the spectrum with a 2'^'* order polynomial and 
divide the spectrum by this fit. This effectively removes 
power on broad-band scales (i.e. ~ 2000A resolution). 
The spectra are then convolved with a 200A top-hat filter 
to make a smoothed version and are then divided by this 
smoothed version. Thus only high-pass spectral informa- 
tion is retained. The x^ is then the fit between the high- 
pass model and the high-pass data, using a suitable noise 
model. This is 'FOM B' in the nomenclature of BG02. 
Another x^ is computed from the low-pass 200A resolu- 
tion spectra (FOM A). Finally we also compute a x^ fit 
between SDSS ugriz colors and PEGASE model colors 
which is illustrated in the lower panel of Figure 4. This 
final FOM is marginalized over a range of extinction corre- 
sponding to about Ay = 0.5-2 mags because of the extinc- 
tion dependency of this FOM. The upper panel of Figure 4 
illustrates the steps in this procedure for data and example 
model spectra (both bad fits and good fits). 

In general, we use only the high-pass FOM or a weighted 
combination of the FOMs with less weight being given to 
the low-pass spectral / broadband color information. Ex- 
amples of residuals for the high-pass FOM are shown in 
Figure 5. The advantage of this approach is that low-pass 
systematic uncertainties such as flux calibration errors and 
dust extinction are effectively exchided. The disadvan- 
tage is that broadband color information is also mostly 
excluded from the analysis. In practice we find this makes 
little difference since most of the information in the col- 
ors is already encoded in the spectral breaks and line in- 
dices. Additionly, we note that regions near strong emis- 
sion lines are excluded from the analysis as we are inter- 
ested in star-formation history, not the instantaneous SFR 
that the emission lines are most sensitive to. 

To estimate the errors, we follow the approach of BG02 
and divide the sample up into 10 contiguous sky areas. We 
then use the variance between cosmic spectra from differ- 
ent sky areas to estimate the errors, this has the advan- 
tage of including systematic effects. Simply taking the raw 
variance will overestimate the true errors, whilst dividing 
it by 10 will not properly take into account systematic ef- 
fects. As a compromise, we estimate the confidence limits 
with Monte Carlo simulations each time drawing 5 random 
entries from the 10 regions. This will in principle still over- 
estimate the errors by ^2 but is robust against systematic 
effects. In addition, we add 2% random errors to the 
models (with power on all scales from high-pass spectra to 
broad-band photometry) . 

3.2. Results 

BG02 found the 2dFGRS cosmic spectra exhibited no 
unique solution but there was a broad curved degener- 
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Fig. 4. — Determination of figures of merit (FOMs). (loft column is a poor fit model, right column is a good fit model). Upper: spectral 
FOMs. The normalized data and and example model spectrum (step 0, offset by +0.5) are adjusted by a low-order polynomial to remove 
power on broad-band scales and some emission lines are cut out to remove nebular components. This forms an adjusted set of spectra 
(step 1). To compare these, the spectra are further divided into low-pass and high-pass components (step 2 A and B, offset by —0.3 and 
— 1.0). The FOMs are determined from the difference between the model and data for each component. Lower: photometric FOM. The 
broad-band fluxes are coadded without fc-corrections to form the data ^{ugriz) magnitudes where z is the mean redshift of the slice. The 
model magnitudes are calculated by integrating the model spectrum through the redshifted bandpasses. Both are normalized by subtracting 
the mean of ugriz. The FOM is formed from the minimum difference after comparing the data with 10 different extinction versions of the 
model magnitudes (dotted lines). The squares represent the fiducial extinction model. 
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Fig. 5. — Example of residuals from the high-pass analysis: residual squared versus wavelength. The middle panel represents a good-fitting 
model while the outliers show models rejected at very high formal confidence (> 99.99%) in comparison with the spectrum from SDSS 
region A. 
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acy surface in the a — (3 plane. However it did allow the 
defining of a common solution where all the different mea- 
surements of a — /3 agreed. As discussed by Hogg (2002) 
the surveys to z = 1 favor /? ~ 3. Direct UV observations 
of high redshift galaxies give a ~ — 1 (Madau et al. 1996), 
i.e. a rapid decline, but once these are extinction corrected 
for dust, estimates range from constant SFR with a = 
(Madau et al. 1996; Steidel et al. 1999) to increasing star 
formation, with a = 1 for 1 < z < 1.5 (Cowie et al. 1999). 
The problem of the effect of dust extinction, and the sam- 
ple bias effects such as luminosity and surface brightness 
selection, is that increasing amounts of UV are missed as 
one goes to higher redshifts. It has recently been claimed 
that very large amounts of UV and hence star formation 
are missed at high redshifts (Lanzetta et al. 2002). So the 
values of a are really lower limits. This is where the cos- 
mic spectrum constraints become the most useful: only a 
narrow range is allowed and high a — /? is excluded. The 
conclusion of BG02 was that there was a concordance SFH 
of 1.8 < /3 < 2.9 and -1.0 < a < 0.7 that agreed with the 
various determinations at the 95% level. 

We compare SDSS and 2dFGRS derived cosmic SFHs by 
plotting both of them in the a — (3 plane (with r = 1.1). 
This is shown in Figure 6 for the two redshift ranges A 
and C. Contours are shown at the 90% confidence level 
with only the high-pass FOM applied. In all cases the 
general shapes of the allowed, degenerate regions are sim- 
ilar, indicating qualitative agreement between the kind of 
SFHs permitted by SDSS and 2dFGRS cosmic spectra. At 
low redshift (region A), where we expect the samples to 
have the least luminosity bias, the agreement is good and 
well within the formal 90% confidence limits. The differ- 
ence amounts to A/? ~ 0.25 for a = and fixed chemical 
evolution. This leads to our first important conclusion: 
the effects of sample selection in the 2dFGRS paper were 
indeed small at low redshift. 

If we look at the higher redshift contours (region C) 
then there is an obvious shift in the contours from their 
low-redshift positions, which is the effect of the aforemen- 
tioned luminosity bias and sample selection. This becomes 
clearer if we can plot the trend with redshift, to do this in 
Figure 7 we fix a = and plot /? versus redshift for 2dF- 
GRS and SDSS redshift slices. The plot shows 95% con- 
fidence regions. Here, we apply a weighted combination 
of low-pass/broadband and high-pass FOMs but qualita- 
tively the results are the same if only the high-pass infor- 
mation is used. At low redshift, the two regions converge 
and overlap considerably around 2 < /? < 4 whilst at high 
redshift the two diverge, which we interpret as the effects 
of differential luminosity bias in the red versus blue selec- 
tion. For the red sample the derived (3 for a slice increases 
considerably approaching /3~6at2;~0.2, indicating 
an increasing dominance of luminous but old red galaxies. 
In the blue sample this is somewhat counteracted as the 
blue selection below 4000A can include more young star- 
forming galaxies at these redshifts leading to lower values 
of/3~3at2;'^0.2. This can be seen directly in the colors: 
if we match up 2dFGRS redshifts with SDSS photometry 
we also see much bluer u — r colors for a ^ > 0.2 and 



6j<19.45 sample (matching the 2dFGRS selection) than 
for a 2; > 0.2 and r < 17.7 SDSS selection. The effect 
is strongest for z > 0.2 where the 6j band is pulled be- 
low the 4000A break. We note that at low-redshift the 
u — r distribution is bimodal with a red-type and blue- 
type branch (Baldry et al. 2003) whereas at z > 0.2 the 
blue peak is gone in the r selected SDSS sample but is 
still there in the 6j selected 2dFGRS sample. The effect of 
producing a more constant f3 is somewhat fortuitous. One 
can try making a blue-selected subset of the SDSS spectra, 
however to be complete for all galaxy types it has to be 
much brighter (approximately 6j<18.1) and the redshift 
distribution is too low to show significant effects (98% of 
galaxies have z < 0.15 with a median redshift of 0.07). 
Howciver the important point is that for the lowest red- 
shift slices shown here we have a good approximation to 
a complete volume limited sample and that the selection 
function becomes small. For a = we find from the joint 
overlap of 2.5 < /? < 3.5 (region A, 95% confidence). If a 
is increased then (3 comes down. 

We have demonstrated by comparison with the com- 
pletely independent SDSS sample that the effect of lumi- 
nosity bias, one of the outstanding issues in the 2dFGRS 
result of BG02, is small. We note the remaining problems 
of blue versus red samples could be improved, particularly 
at high redshift, by a multiwavelength selected sample, 
for example all galaxies jointly in each of ugriz down to 
a magnitude limit in each band which matched the esti- 
mated cosmic spectrum at some redshift. 

The remaining important bias issue to consider is that of 
aperture effects. The SDSS sample allows us to investigate 
this for the first time using the colors. In Figure 8 we plot 
the difference in the average u — r color between 'model' 
SDSS magnitudes and 'aperture' SDSS magnitudes which 
we have computed by integrating the database image pro- 
file information to a 3 arcsec diameter aperture. This ap- 
proximates the SDSS 'fiber' magnitudes but has the ad- 
vantage of being PSF corrected. We choose m — r as it 
straddles the 4000A break and should be a color that is 
very sensitive to changes in stellar populations (Strateva 
et al. 2001). The 3 arcsec diameter aperture matches the 
spectroscopic aperture. The model magnitudes are de- 
termined by using the best fit of de Vaucouleurs and ex- 
ponential spatial profiles to calculate a 'total' magnitude. 
We do not use the SDSS Petrosian magnitudes, although 
these; work well for total magnitudes in griz bands we have 
determined there are significant problems in the use of u 
Petrosian magnitudes. 

It can be seen from Figure 8 that there is no significant 
difference in the average u — r between model and fiber 
and no trend with redshift for z > 0.05. (We note however 
that at high redshift 2; > 0.1 the SDSS becomes dominated 
by early-type galaxies for which we naturally expect small 
aperture effects in the colors). For z < 0.05, there is a 
small aperture effect amounting to A{u — r) < 0.05 (re- 
gion A). We stress that for individual galaxies there is a 
large dispersion in this statistic because aperture bias is 
a significant problem for individual galaxies. Only in the 



Specifically, the fiux can be summed over large Petrosian radii causing significant magnitude errors because the u-band fiux is so weak and 
objects are so close to the background level in the firames. This can be seen for example by comparing the distribution oi u — r with camera 
column. This is a problem with u magnitudes from version 5.2 of the SDSS PHOTO software and may be improved in the next version. Please 
contact the authors for further information. 
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Fig. 6. — Comparison of star-formation histories derived from the SDSS and 2dFGRS samples in the a — f} plane for cosmic spectra from 
redshift regions A and C. The contours represent 90% confidence in our analysis using only the normalized high-pass spectra. For a given 
a = 0, the difference between the SDSS and 2dFGRS results amounts to A/3 ~ 0.25 for the low-redshift samples and A/3 ^ 1 for the high 
redshift samples. 
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Fig. 7. — Comparison of star-formation histories derived from different redshift slices. For each slice we plot the 95% confidence region 
on /3, for a = 0. The upper region represents the SDSS contours (solid bounding lines) and the lower region represents the 2dFGRS con- 
tours (dashed bounding lines). The variation in this plot docs NOT represent any cosmic evolution, rather it shows primarily the effect of 
luminosity-selection bias in measuring cosmic SFH (represented by /3). The divergence is due to the red versus blue selection in the samples 
as one goes out to high redshift. At low redshift where most of the cosmic light is being counted by both surveys the values of recovered /3 
converge. The A,B,C lines show the redshift ranges defined in Table 1. Note that the results here differ slightly from Fig. 6 because here we 
have marginalized over metallicity in the range 0.5-1.5 Zq. 
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Fig. 8. — Plot of (u — r)APER — {u — r) model average colors for a fixed luminosity range in the SDSS sample versus redshift to illustrate 
aperture effects. (See text for definition of our aperture magnitudes). The x symbols represent the median, the + symbols the mean and the 
vertical lines denote the standard deviation (in redshift bins of 0.005). (We plot the latter instead of the standard error on the mean in order 
to illustrate the scatter). The median and mean values show a small offset at z < 0.05 but it drops to zero at higher redshift. The dispersion 
remains large indicating for individual galaxies, as opposed to the mean, aperture effects remain a problem. We also note that the sample is 
not representative for z > 0.1 (see text). 



mean does the aperture bias appear small. Indeed it ap- 
pears aperture bias from using this low-rcdshift sample is 
smaller than the luminosity bias from using higher-redshift 
samples (region B or C). Given there is no significant aper- 
ture effect in the colors, how docs this translate into stel- 
lar populations? As we know color is degenerate with age 
and metallicity, so either the stellar populations are the 
same or there is some fortuitous cancellation in age and 
metallicity gradients. The latter seems unlikely. The best 
evidence is that aperture effects are not significant for the 
cosmic spectrum, i.e. on average. To investigate this in 
more detail would require two dimensional resolved spec- 
troscopy of a large number of galaxies. We do note that 
for individual galaxies, for example a face on spiral with a 
significant old bulge population, aperture effects could be 
significant. However, there are many types and different 
orientations of galaxies in a magnitude- limited survey and 
the average fiber placement on an average galaxy is more 
robust to aperture effects. 

The final value of /3 is still in the range of previous 
estimates. From SDSS alone: for a = 1 we constrain 
2 /3 ^ 3.5, and; for a = we constrain 2.5 < /3 < 4. 
These ranges cover metallicity degeneracy from about so- 
lar down to half-solar with approximately 95% confidence 
limits. The /3 values are close to the 2dFGRS results of 
BG02. The low-redshift values from the luminosity den- 
sity measures in UV, nebular lines, far-infrared and radio 
compiled by Hogg (2002) give a mean /3 = 3.1 ± 0.7 (Icr), 
with a spread in most of the measurements (5 out of 9) 
of 2.9 < /3 < 4.5. The cosmic spectrum and luminosity 



density measurements are entirely consistent, within their 
own intrinsic uncertainties, and suggest that /3 ~ 2-4. 

3.3. Tests of high-redshift star formation 

As well as the a — jS parameterization we can test more 
complex models of SFH. Recently (Lanzetta et al. 2002) 
proposed that direct high-2; surveys were missing a large 
component of high redshift light based on a comparison 
with the specific SFR intensity distribution and the col- 
umn density of Lyman-a absorbers which they claimed 
were related. This is an ideal situation to test with the 
cosmic spectrum because even if objects are obscured at 
high redshift the stars produced must end up in objects 
in the low-redshift SDSS data, unless of course they form 
some new local population yet to be detected. We note 
that Lanzetta et al. proposed these were the ancestors of 
todays galaxies. They proposed three possible models us- 
ing different missing-light corrections, from different meth- 
ods of estimating the SFR intensity distribution function, 
which we term 'LOW, 'MEDIUM' and 'HIGH' based on 
their z > 2 SFRs (Figure 9). All have a low-z slope equiv- 
alent to /3 « 1.2. Note we do not use the a — /3 parameteri- 
zation, rather we fit the SFHs given by Lanzetta's figure 4 
directly with a break at z « 2 (see Figure 9 where we re- 
plot the Lanzetta figure to illustrate the time dependence 
and our fit to Lanzetta's points). The results of fitting 
these SFHs are given in Table 2. We use a weighted com- 
bination of FOM A with FOM B in the ratio of l-to-5 
based on uncertainties determined from the Monte-Carlo 
error estimation from our different sky regions. For refer- 
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Fic;. 9. Three possible cosmic SFHs from the analysis of Lanzetta ct al. (2002). The symbols represent their data while the solid lines 
are our fitted to their SFH that we use in our test. The variation is principally at very early times, as illustrated in the right panel. The 
timescale shown is for the standard A cosmology with Hq = 70. Note that we only use the average values for each Lanzetta et al. model, i.e. 
we do not consider their error bars. 



cncc wc also give the significance levels of a fiducial a = 0, 
/3 = 3 model and a a = 1, /? = 2 model. The first gives 
a best-fit model with about half-solar metallicity and the 
second with about solar metallicity. 

The HIGH model is rejected at > 99% confidence for 
cosmic spectra from both the SDSS and 2dFGRS surveys. 
In fact, the HIGH model only fits at this level if the av- 
erage metallicity becomes rather low {Z < 0.5Zq), if we 
restrict the range to > 0.5.^0, the model is rejected very 
severely. The LOW model is also rejected by the SDSS 
data because of the low /?, though more marginally. None 
of these models fit the cosmic spectra as well as our fidu- 
cial a = 0, P = 3 model which is well motivated by the 
high redshift luminosity density measurements. Thus it 
appears the cosmic spectrum does not provide any strong 
evidence for largo amounts of missing light in the current 
high redshift census. Quantitatively, no more than about 
85% of stars formed at ^; > 1. This is slightly higher 
than the BG02 upper limit of 80%. The differences are 
that wc include the SDSS data, include only the 2dFGRS 
low-rcdshift range, give some weight to the low-pass com- 
ponents of the cosmic spectra and allow ^form = 10 for the 
Lanzetta models. It is of course possible that some of our 
assumptions may be wrong, such as a universal IMF; how- 
ever the data can be well fitted by models which assume 
a Universal IMF. The slope of the high-mass IMF can be 
constrained by cosmic spectra, if near-IR data is included. 
This will be addressed by a forthcoming paper (Baldry & 
Glazebrook 2003). 

4. THE ABSOLUTE COSMIC SPECTRUM 

We have shown that we have derived a normalized cos- 
mic spectrum which is consistent between SDSS and 2dF- 
GRS surveys, is close to being volume limited and is robust 



against aperture effects (by comparing with SDSS u — r 
colors). An advantage of the SDSS survey is the excel- 
lent spectrophotometry as standard stars are included on 
each plate. In particular we would expect that the rela- 
tive fluxing with wavelength of the SDSS cosmic spectrum 
should be much more accurate than for 2dFGRS (BG02 
indicate 5-10% errors are likely in the latter). We do in- 
deed find that the SDSS cosmic spectrum is in agreement 
with the 2dFGRS cosmic spectrum once a good jit spec- 
trophotometric correction has been applied to the latter 
(sec BG02 for details of this procedure, essentially it uses 
a good fit high-pass model to set the fluxing). 

So far we have been dealing with the normalized cosmic 
spectrum, and principally fltting to the high-pass spectral 
information. However we can put the cosmic spectrum 
on an absolute luminosity scale by normalizing to the r- 
band luminosity density. For consistency, we calculate this 
ourselves from a more recent large-scale structure (LSS) 
sample by calculating the r luminosity function using 
Kiurvcy/Knax Weighting and fc-correcting to the rest-frame 
r-band (kcorrect vl_ll, Blanton et al. 2002b). The LSS 
sample is similar to our cosmic-spectra sample but has a 
well defined area, includes nearest-neighbor redshifts to re- 
place galaxies missed due to fiber collisions and has stricter 
limits on the selection magnitude of 14.5 < r < 17.5 (taken 
from scunplelO described by Blanton et al. 2002c). A 
Schechter function fits very well and the total luminosity 
density is calculated from the analytic total integration 
(Figure 10). We note there is a small excess of galax- 
ies above the Schechter function flt at low luminosities 
Mr > —18.5, a disrepancy that has been noted by other 
authors (Madgwick et al. 2002) but this has negligible ef- 
fect on the flnal luminiosity density (see the lower panel 
of Figure 10). Regions A and B sample better the faint 
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Table 2 

Comparison of SFH models with cosmic spectra 



SFH Model 


SDSS confidence" 


2dFGRS confidence"* 


a = 0, 13 = 3 


<68.3 


68.3 


a = 1,13 = 2 


80 


<68.3 


Lanzetta LOW 


98 


68.3 


Lanzetta MED 


90 


90 


Lanzetta HIGH 


99 


99 


Lanzetta HIGH Z > 0.5Zq 


99.99 


99.99 



° Percent rejection for models marginalized over metallicity. The 
models assume a universal Salpeter IMF and are compared with 
the low-redshift range cosmic spectra from the surveys. Lanzetta 
models are integrated for z < 10. 



and bright end of the luminosity function respectively, the 
luminosity density in these two redshift ranges agrees to 
±6%. For our final r luminosity density we take the value 
calculated from regions A+B (0.015 < z < 0.08), this is 
1.94±0.20x 10^ h Lq Mpc-3 (j + 2.51og/i = -16.1±0.1 
mags where j is the absolute magnitude of the integrated 
light per Mpc^). We estimate that the systematic uncer- 
tainties (redshift ranges, Schechter fitting) are about 5- 
10% and quote a 10% error bar. The final value agrees 
well with the more sophisticated luminosity density calcu- 
lation of Blanton et al. (2002c). 

We normalize our cosmic spectrum for region A (where 
it is the least affected by luminosity bias) to this r luminos- 
ity density by integrating the flux of our cosmic spectrum 
through the r-band filter profile (Gunn et al. 1998) and 
calculating a scaling factor. We note that this will correct 
to first order for any residual luminosity bias, as the lu- 
minosity density is calculated by integrating the Schechter 
function fit to zero luminosity. 

This absolute cosmic spectrum (units: Watts 
Mpc~^) is shown in Figure 11 (smoothed slightly for the 
plot) and tabulated in Tables Al and A2. We can check 
the reliability of this by comparing the spectrum with lu- 
minosity densities computed in other broad bands, these 
should correspond to smoothed estimates of the cosmic 
spectrum. In the figure, we overlay the SDSS ugriz lu- 
minosity densities from Blanton et al. (2002c) the 6j 
luminosity density from Norberg et al. (2002), and the 
ugriz luminosity densities calculated by using the cosmic- 
spectrum sample directly (i.e. combining the photometric 
ugriz fluxes with the same weighting as the spectra). 

There is excellent agreement within the errors. We note 
that the luminosity densities given by Blanton et al. (2001) 
are 20-30% discrepant in u and g, but the more recent de- 
termination shown here from the Blanton et al. (2002c) 
sample of 150,000 galaxies are in much better agreement 
with this work, with Cross & Driver (2002), Norberg et al. 
(2002) and the revised SDSS work of Yasuda et al. (2001). 
We note some minor differences in analysis: Blanton et 
al. uses a maximum likelihood method which allows for 
evolutionary effects (in a set of magnitude-limited sam- 
ples in each band separately) whereas the Cross & Driver 



analysis used the l/Vmax method, and the redshift ranges 
did not quite match (although both were more or less 
< z < 0.15). 

Finally we use our physical cosmic spectrum to derive 
some interesting quantities. The mass/light ratio of our 
best-fit models is 3.7-7.5 Mq/L© in the r-band (assuming 
Av « 1). The range corresponds to the range from a low- 
mctallicity a. = model to the solar metallicity MEDIUM 
Lanzetta model. Thus from the r-band luminosity den- 
sity we can derive the cosmological mass density in stars: 
^stars/i = 0.0025-0.0055. Of course this is highly depen- 
dent on the assumed Salpeter IMF as the galaxian light is 
dominated by the most luminous stars. For example the 
mass-to-light ratio for an IMF with slopes (—0.5,-1.35) 
broken at 1 Mq is about 60% of the unbroken Salpeter 
IMF (—1.35). We also computed ilstars^ for the Kennicutt 
(1983) IMF (-0.4, -1.5), redoing the fitting, and obtained 
0.001 0.002. 

Similarly our best fitting range of models also gives us a 
mean cosmological SFR today, the value range is 0.01-0.04 
h Mq yr~^ Mpc^''. For a particiflar choice of a. and /? the 
range is narrowed as there is only a remaining degeneracy 
with metallicity. For example if motivated by the high- 
redshift luminosity density measures we choose a = 0, 
/3 = 3 we find the SFR today is 0.025-0.03 h Mq yr'^ 
Mpc~^. It is also interesting to look at the emission lines 
in the spectrum, by subtracting the stellar 'continuum'^^ 
and integrating the fiux we can derive the local luminosity 
densities of a whole set of useful lines: [Oil], H/3, [OIII], 
Ha, [Nil] and others. We can do a continuum subtraction 
by subtracting our best fit SFH stellar model, although 
the resolution is different the subtraction is still close to 
zero except near the nebular lines. We used the best fit 
high-pass model, so as to get the best match to absorption 
features, and normalize it to the physical cosmic spectrum 
using a smoothing filter. Wc show the example contin- 
uum subtracted spectrum in Figure 11. The subtraction 
is excellent in the 4000A-8000A region, the principal dis- 
crepancy is for wavelengths > 8OOOA where the Calcium 
triplet absorption lines are not well fit by any of our mod- 
els; this effect can also be seen in Figure 5. 

The line luminosity densities are simply calculated by 



This is based on a sample of 80 000 galaxies and represents an update on Blanton et al. (2001) 
By 'continuum', we mean stellar continuum and absorption features. 
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M, (h = 0.7) 

Fig. 10. — Luminosity function determined from redshift regions A+B. Upper: The number of galaxies Mpc^'^mag^^ versus magnitude 
with a Schcchtcr function fit (solid line). The vertical bars represent ±3(t Poisson errors. Lower: A linear scaling of luminosity (lO^'^'^/-^) 
versus magnitude, where Ai is the total absolute magnitude Mpc~'^ mag~^. This shows that the integration of the Schechter function is a 
good approximation to the luminosity density. This plot uses Petrosian magnitudes but the result is similaj: if Sersic magnitudes (Blanton 
et al. 2002c) are used; j +2.51og/i = -16.11. 



integrating tlio cosmic spectrum for each hne from A — 8A 
to A + 8 A. Tlio box width is chosen to be 3 x the typical 
line FWHM so the fluxes are very close to total. 

The resulting line luminosity densities are given in Ta- 
bic 3. Wc quote values for region A and B, as the latter 
includes [Oil] and for comparative purposes. Both are 
normalized to the same r luminosity density. There are 
two sources of error: firstly the imperfection of the contin- 
uum subtraction. Visually we checked the region around 
each line, particularly the Balmer lines, and found no ev- 
idence of significant residuals, i.e. the plot appears as an 
emission line with zero continuum plus 'noise' due to the 
mismatch in resolution between raw data and model spec- 
trum. To quantify this for each line we similarly extract 
two regions, either side of the line, of the same width where 
there is only blank continuum, i.e. typically a 40-60A off- 
set. (These are optimized for each line particularly in the 
crowded Ha/ [Nil] region). These give typically the same 
error and we quote the maximum in the table. There is 
also an additional systematic error due to the uncertainty 
in the r luminosity density which we estimate as 10%. 

The line luminosity densities lead to some interesting 
results: first the Balmer decrement Ha/H/? is 6.4±0.9 (re- 
gion A but region B is similar), if we take an unreddened 
case B recombination value of 2.86 (Hummer & Storey 
1987) and a Milky Way dust law from Pei (1992) (the 
SMC law gives similar numbers as the two are close in the 
optical) we derive a nebular extinction Ay = 2.4 ± 0.4. 
The Ha luminosity density thus requires a correction fac- 
tor of x(5.8 ± 1.8). This is consistent with other workers 
findings — for example we have re- analyzed the sample of 
Gallego et al. (1995) looking at the Ha correction factor 
as a function of luminosity and find a volume averaged 



value of x4.6. Thus it makes little difference whether one 
dereddens before taking the mean or simply dereddens the 
mean (which is the approach we are taking here with the 
cosmic spectum) even for these large nebular extinction 

values. 

The dereddened Ha luminosity density is thus 4.1 ± 
1.3 X 10^^ W Mpc~^ for h = 0.7. Converting to cgs units 

and h = 0.5 this gives 2.9 ± 0.9 x 10^^ ergs s'^ Mpc'^ 
which is the same as that found by from the Canada France 
Redshift Survey (Tresse & Maddox 1998b) at low redshift 
{z 0.2). It is twice as high as that found by the objective 
prism survey of Gallego et al. (1995). 

From our models we can also work out the conversion of 
Ha luminosity into SFR (from the number of ionizing pho- 
tons), the mean conversion factor is 1.36 x 10^^ ergs s"""^ 
Mq"^ yr and the range in the models is ±15% (it varies 
with metallicity). This allows us to derive a SFR density 
today of 0.030 0.056 h Mq yr^^ Mpc"^ It is remarkable 
to note that this is entirely consistent with the derived in- 
dependently from the star-formation history fitting which 
illustrates the validity of our general approach. The Ha 
line independently measures the instantaneous SFR from 
young ionizing stars, whereas the SFH fits a model to the 
stellar continuum and is averaged over several Gyr. It is 
interesting to note that given the cosmic spectrum to get 
a SFR from SFH fitting much lower than 0.03 h Mq yr~^ 
Mpc^'' requires high /3, high a AND a low metallicity to- 
day. 

We note that the Gallego et al. (1995) sample also gives 
a lower SFR density than the UV sample of Treyer et al. 
(1998). Both Treyer et al. (1998) and Tresse & Maddox 
(1998b) are measured at 2 ~ 0.2, whereas the current sam- 
ple is 2; 0.05 which is comparable to the redshifts of the 
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Fig. 11. — Absolute cosmic spectrum from SDSS region A, i.e. luminosity density per unit wavelength in physical units (dotted line, 
h = 0.7). This has been resampled at 10 Angstroms for display clarity. The normalized cosmic spectrum (Fig. 3) was scaled to match the 
r-band luminosity density (Fig. 10) shown by the triangle. The squares show the relative luminosity densities, i.e. flux-weighted observed 
colors, from region A scaled to the same absolute luminosity density (note they are offset in wavelength to account for the mean redshift of 
0.035). The other points represent different estimates of broad-band luminosity densities: asterisks show the redshifted (ugriz) filter results 
of Blanton et al. (2002c) and the diamond shows the f)j point of Norborg et al. (2002). The lower spectrum shows the stellar-component- 
subtracted cosmic spectrum described in the text (shown at the same resolution as above but line luminosities were estimated from the full 
resolution spectrum as in Fig. 3). This is still in physical units for the line fluxes, i.e. the ordinate units still applies. 



Gallego et al. (1995) sample. It is only the latter which is 
discrepant. The most likely explanation is that the Gal- 
lego et al. (1995) survey covers only 1/5"^ our sky area and 
is only sensitive to star-forming galaxies as it is emission 
line selected. Our survey represents a 'cosmic average' and 
galaxies contribute to the final cosmic spectrum even if the 
Ha flux is not detectable in individual galaxies. We note 
that Ha is only sensitive to transient star-formation over 
~ 20 Myr whereas a UV sample like those of Treyer et al. 
(1998) effectively average over longer times scales of up to 
~ 1 Gyr (see Glazebrook et al. (1999) for a discussion of 
this). 

A final comment on the other line ratios: the observed 
Ha/[OII] line ratio of 2.1 (region B) is entirely consistent 
with the median found by the sample of galaxies observed 
by Kennicutt (1992). The other line ratios are also entirely 
consistent with those from star-forming galaxies (Veilleux 
& Osterbrock 1987) including the weak [01] line, indicat- 
ing that the AGN contribution to the cosmic spectrum is 
indeed negligible. It can at most be only a few percent 
according to the models of Kewley et al. (2001). Most of 



the optical light of the Universe does indeed come from 
stellar nucleosynthesis. 

5. SUMMARY AND CONCLUSIONS 

Computing the light in the Universe today from the 
SDSS and 2dFGRS surveys allows us to derive a cosmic op- 
tical spectrum, determine its robustness, and make more 
accurate determinations of allowable star-formation histo- 
ries of the Universe. In particular: 

1. We find a range of solutions 2 < /3 < 3 and 
< a < 1 which are consistent with the cosmic 
spectrum and determinations of the evolution of the 
luminosity density in various bands. 

2. 'Fossil Cosmology' and direct cosmology agree: i.e. 
the SFH inferred from the local Universe agrees 
with that measured by luminous emission at high 
redshift. Again the Copernican Principle, that we 
are in no special place in the universe, is demon- 
strated by observational astronomy. 
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Table 3 

Luminosity density in various lines derived from the cosmic spectrum 



Line 


Luminosity Density (10^" W Mpc 
Region A Region B 


[Oil] 


3727A 




28 ±4 


H/3 


486lA 


11±2 


9± 1 


[OIII] 


4959A 


6±2 


3±2 


[OIII] 


5007A 


22 ±2 


13 ±2 


[01] 


6300A 


2± 1 


2± 1 


[Nil] 


6548A 


11± 1 


11± 1 


Ha 


6563A 


70 ± 1 


59 ± 1 


[Nil] 


6583A 


30 ± 1 


29 ± 1 


[SII] 


6716A 


13 ± 1 


11± 1 


[SII] 


6731A 


10 ±2 


8±2 



Notes: (i) Values are quoted for h — 0.7. Errors on each value are 
for the continuum subtraction error; (ii) Additional errors of ±10% 
should be added to allow for the systematic uncertainty in the r 
luminosity density. 



3. There is good agreement between SDSS and 2dF- 

GRS derived cosmic spectra and SFHs at low- 
redshift, where we expect luminosity biases to be 
minimal, despite the difference in aperture. We 
conclude the result is robust against luminosity se- 
lection effects. 

4. The excellent photometry of SDSS allows us to test 
for aperture effects; for average quantities such as 
the mean color and hence the cosmic spectrum we 
find the aperture effects are not significant (though 
they are still important for individual galaxies). 

5. Due to the excellent spcctrophotomctric quality of 
the SDSS data we can make an absolutely cali- 
brated cosmic spectrum for a close to volume lim- 
ited sample. 

6. None of the SFH scenarios proposed by Lanzetta 
et al. (2002) fit the cosmic spectrum as well as more 
standard models. There appears to be no com- 
pelling evidence for missing star-formation at high 
redshift from the cosmic spectrum. 

7. Typically good fits to the cosmic spectrum (using 
consistent chemical evolution) give a final metallic- 
ity of 0.5-1 Zq. 

8. We find the stellar population of the universe today 
has a best fit r-band mass/light ratio (for Salpeter 
IMF) of 3.7 7.5 M0/L0, which given the r-band 
luminosity density of j + 2.5 log /i = —16.1 gives 
^^stars/t = 0.0025-0.0055 (and a factor of 2 lower for 
a Kennicutt IMF). 

9. By fitting the best stellar population model and 
subtracting, we can derive a whole set of nebular 
line luminosity densities for our cosmological vol- 
ume. In particular we find the local dereddened 
Ha luminosity density is twice as high as found by 



Gallego et al. (1995) and similar to that found by 
Tresse & Maddox (1998b). 

10. We find the SFR of the Universe today (for Salpeter 
IMF) is 0.03-0.04 h Mq yr"! Mpc'^, and agrees be- 
tween models fit to the stellar population and that 
derived from the Ha luminosity density of the same 
sample. 

5.1. Future Work 

One avenue which it is clearly possible to explore is how 
the SFH varies with the luminosity of galaxies. The cos- 
mic spectrum represents a binned approach, i.e. looking at 
the ensemble stellar population in a volume and inferring 
the SFH. We can not specify which galaxies the stars were 
in at earlier times, i.e. it is not sensitive to merging. We 
can extend this approach by computing the cosmic spec- 
trum in luminosity bins, i.e. estimating the SFH of stellar 
populations as a function, approximately, of the mass of 
the galaxy they end up in today. Does the SFH depend 
significantly on this? Given known color-luminosity rela- 
tionships which extend across a range of Hubble types (e.g. 
Gavazzi 1993, figure 4) obviously we expect it too. 

In fact we have seen already some evidence for this in 
the current paper, where the selection effect of the magni- 
tude limit causes the high redshift part of the SDSS survey 
to have higher luminosities and steeper values of /3. Lumi- 
nous galaxies on average contain older stellar populations. 
In the next paper in this series (Baldry et al. 2003) we 
present a detailed analysis of this differential SFH in to- 
day's Universe. 

One limitation of this work is the assumption of a uni- 
versal Salpeter IMF. The high-mass IMF can be con- 
strained by cosmic spectra, if near-IR data is included. 
This will be addressed by a forthcoming paper (Baldry et 
al., 2003). 

Another limitation of the current work is the depen- 
dence on relatively low resolution (20A) models of evolu- 
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tionary synthesis, whereas the SDSS data has 2-5 A res- 
olution. The SDSS cosmic spectrum resolves many low 
equivalent width metal lines which are not exploited by 
the current low resolution analysis and will help to break 
the agc-mctallicity degeneracy. We arc working toward 
being able to construct high-resolution models to resolve 
further the question of the star-formation history of the 
Universe. 
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Rest Wavelength 

A 


Luminosity 
103" W Mpc"3 


Continuum subtracted luminosity 
lO^o W Mpc~^ 


6556.6 


a A cm n 


U. oU (DU 


6558.1 


7.10810 


96345 


6559.6 


9.73218 


3.62092 


6561.1 


15.42100 


9.34124 


6562.7 


19.65739 


13.61270 


6564.2 


16.94022 


10.92763 


6565.7 


10.88128 


4.90330 


6567.2 


7.43149 


1.47872 


6568.7 


6.54169 


0.62419 


6570.2 


6.34205 


0.44839 



Note: Only a sample 10 rows are shown in a wavelength region near Ha. The full table covering 3764A-8697A is available 
from the electronic edition. 



Table A2. Cosmic Spectrum for Region B. 



Rest Wavelength 

A 


Luminosity 
1030 W A-i Mpc-3 


Continuum subtracted luminosity 
1030 w A-i Mpc-3 


3723.7 


4.10668 


1.17038 


3724.6 


5.54068 


2.62494 


3725.4 


6.87602 


3.97530 


3726.3 


7.58438 


4.70544 


3727.2 


8.05617 


5.19847 


3728.0 


8.50139 


5.66404 


3728.9 


8.14010 


5.32336 


3729.7 


6.53688 


3.74197 


3730.6 


4.65805 


1.85869 


3731.4 


3.45462 


0.63813 



Note: Only a sample 10 rows are shown in a wavelength region near [Oil]. The full table covering 3636A-8458A is available 
from the electronic edition. 



